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Coulomb problem for a Dirac particle in flat Minkowski space 
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Abstract 

It is shown that there exist several ways to treat the quantum-mechanical Coulomb 
problem for a Dirac particle in flat Minkowski space with the help of the Heun differential 
equation, Fuchs's equation with four singular points. When extending the problem to curved 
spaces of constant curvature, Lobachevsky H 3 and Rimann S3, there arise 2-nd order differ- 
ential equations of the Fuchs type with 6 singular points, a method to get relevant equations 
with five singular points has been elaborated. 

Let us start with a Dirac equation written in spherical tetrad of flat Minkowski space (more 
detail see in pQ) 



i 7 dt + i (7 3 d r + 



7 1 a 31 + 7 2 a 32 , 



+ - S 6 

r 



in 



= , 



ids + i c 12 



sml 



In that basis, spherical solutions are constructed on the base of a substitution (Wigner's D 
functions are noted as D J _ m CT (0, 9, 0) = D a ): 
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ejm 



With the use of recurrent relations [2] 

de D +1 / 2 
—m — 1/2 cos 8 



fi(r) D_ 1/2 
/2O") -D+1/2 
Mr) D_ l/2 
/ 4 (r) D +1/2 



a D_i/ 2 - b D +3/2 , 
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b D -3/2 ~ a D +l/2 , 



where 



-to + 1/2 cos 9 
sin 9 

3 + 1/2 
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-1/2 



-bD_ 



-3/2 



aD 



+1/2 , 



we derive (below v = j + 1/2 ) 



-iet 



- U(r) D_ 1/2 
+ f 3 (r) D +1/2 

+ /a(r) £U/2 

- /i(r) D +1/2 



and after simple calculations obtain radial equations 



e/3 - jj-/3 - «-/4 - m/i = , e/ 4 + i^-fi + «-/3 - m/ 2 = 
ar r ar r 



e/i + jj-Zi + «-/2 - rnf 3 = , e/ 2 - £-^-/ 2 - i-fi - ra/4 = 
ar r ar r 

Usual P-reflection operator [3j in Cartesian basis lie = ® -P 



n c . = 
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P, P(9,cf>) = (ir-9, + vr) 



after translation to spherical basis 

n sph = 5(^,0) n c s , - 1 (f > . 

assumes the form 



n 



sp/i 



From eigenvalue equation 



H-sph^jm = n 
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m,a ( 


<f>,9 


0) = 



p . 



,0) = i-iy & (<t>,9,o) 



it follows 



that is 



n = 6 (-l) i+1 , 5 = ±1 , U = S h, h= 5f 2 ; 



-ict 



Allowing for (2a), we simplify eqs. (lb): 



fi(r) £-1/2 
h{r) D +1/2 
& hir) D_ 1/2 
5h(r) D +1/2 



(2a) 



i-r + - ) f + (e + Sm) g 
ar r 



0, 



d v 

h ) 9 

ar r 

where instead of f\ f 2 new functions are used 

/1 + /2 



/ 



To simplify the wave functions 



( e _ 5 m) f = , 



/1-/2 



(26) 



V2 ' 



ejmS 



iV2 ' 

and one can use so-called Dirac (or Johnson- 



Lippmann) operator [3] K (see also [1]), 

£ = - 7 ° 7 3 = 7 °7 3 
Indeed, from the eigenvalue equation K ty e j m (x) = K 9 e j m we get 



7 1 (do + 1/2) + ^— d 4 



K = -5 (j + 1/2) , 5 = ±1 , 

U = 8 fx , f 3 = 5 f 2 , 

which coincides with the above restrictions. 

Transition to a Coulomb problem for a Dirac particle is reached by formal chane the radial 
system (2b) 

, d v . „ ,„ e 

(:r + -)/ + (^ + - + m )3 = o, 



-) 9 



(E + - - m) / = 0. (3) 

Let us perform a linear transformation over functions f(r) and g(r) (its coefficients may 
depend on the radial variable, let its determinant obey an identity a(r)b(r) — c(r)d(r) = 1) 



f(r) = a F(r) + c G(r) , g(r) = d F(r) + 6 G(r) , 

F(r) = b f(r) - c g(r) , G(r) = -d f(r) + a g(r) . (4a) 

Let us combine eqs. (3) as follows: the first equation is multiplied by +6, the second by — c, 
and add the results; analogously, the firs equation multiply by — d , the second by +a, and sum 
the results. Thus we arrive at 

-b'a + c'd + - (ba + cd) + (E + - + m) bd + (E + - - m) ca 
dr r r r 



b'c - be' 26c - (E H h m)6 2 - (E -\ m)c 2 

r r r 



G 
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4- + d'c- a'b - - (dc + ab) - (E + - +m)bd- (E + - - m) ca 
ar r r r 



d'a + da' + -2a<i + (£+- + m)d 2 + (E + - - m)a 



F . 



(46) 



For simplicity, let the transformation (2) does not depend on r, and let it be orthogonal: 

S 
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cos .4/2 sin A/2 
-sin A/2 cos A/2 



which simplifies eqs. (3) 



— H — cos A — m sin4 

dr r 



v e \ 
— sin .4 E — m cos A) G , 

r r J 



cos A + m sin G = ^— - sin A + - + E — m cos A^j F . 



d v 
dr r 

There exist four possibilities (in fact, only two of them are different) 



1) 



1') 



sini+- = 

r r 



sin A = — 



cos 



A = yjl-e 2 /v 2 , 



A 

cos — 
2 



v + \Jv 2 - e 2 
2^ 



sin — 
2 



— sin A = 



sin A = — 



cos 



A 

cos — 
2 



— ^ v 2 — e 2 
2^ 



sin — 
2 



V — \J V 2 


-e 2 


2v 


A=y/l 


— e 2 /v 2 


V + V^ 2 


-e 2 



2v 



(5) 



(6a) 



2) 



2') 



E — m cos A = , cos^l = +— -, sin^l = y 7 ! — E 2 /m 2 , 
^ Im + E A Im-E 



cos ■ 



m + E 
2m ' 



8111 ■ 



E 



2m ' 



-E 1 - mcos^l = , cosA = , sin A = y 7 ! - E 2 /m 2 , 



cos — 
2 



m — E 
2m 



.A m + E 
2 V 2m 



First, consider the case 1). Eqs. (5) takes the form 



(66) 
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d v 



2e 



— H — cos A — m sin A ) F = E - m cos A) G , 

dr r J \ r J 

— — cos A + m sin A ) G = (E - m cos A) E . 
dr r J 

After excluding the function F, we get a second order equation for G 



d v 



dr r 



d v 



— H — cos A — m sin A) cos A + m sin A ] G 



dr r 



2e 



or in more detailed form 



(E - m cos A) ( E - m cos A ) G , 



d 2 <> o 1/ cos A - v 2 cos 2 ^4 
— 7 + E 2 - m 2 + , 



<ir 2 



2eE — 2em cos A + 2mv sin A cos A N . , , 
H | G = 



Taking in mind identity sin A = e/v, thus equation reduces to 



dr 2 



+ E 2 -m 2 + 



v cos A — v 2 cos 2 A 2eE 

n + 



G = 0. 



After changing the variable, x = 2\/m 2 — E 2 r, it reads 



d 2 G ( 1 vcosA{ycosA-\) 
dx 2 



+ 



eE 



4 x' 2 ' yjm 2 - E 2 x 

With the use of a substitution G(x) = x a e bx ip(x), for ip we get 

d 2 (p , , , do? 

x-^- + 2a + 26x 
aar ax 



G = 0. 



+ 



,, 9 1 N a 2 — a — z^cos A (z/cos ^4 — 1) , eE 

6 2 - -) x + i + 2ab + ; 

V 4 ; x Vm 2 - E 2 



When 



+v cos A = \J v 2 — e 2 , 



2 ' 



this equation for 99 becomes simpler 



which is a confluent hypergeometric equation 



eE 



Vm 2 — E 2 



<p = 0, 



To polynomials there correspond the know restriction a = — n, n = 0, 1,2, which gives the 
known energy quantization rule 



eE 



a — 



= —n 



E = 



m 



^ m2 ~ E2 " yfl + e 2 /(n + v 7 ^^) 2 ' 

In turn, from (7) it follows a second order equation for F 



(8c) 



d v . 
cos A + m smA 

dr r 



d v 

. — H — cos ^4 — m sin ^4 ) F 

2e + (E + m cos A)r \dr r 

(-E + m cos A)F , 



(9a) 



or 



, d v . d v . 

(- cos ^4 + m sin AH — H — cos A — m sin A) 

dr r dr r 



d i 
+ — In 

dr 



, d v . . 
,, , , — + - cosA-m sini 

2e + (i? + m cos A)ry ar r 

2e 

( h E 1 + m cos A) (-E + mcosA) F , 



and further 



+ 



H — 2 cos A 



dr 2 

m sin A 

r 

+E 2 - m 2 



E + m cos A 



r 2e + (E + m cos A)r J dr 



d 



2e 



r r[2e + (E + m cos A)r] / r 



(— cos A — m sin A) 



v cos A e 2 — v 2 2eE 



+ 



+ 



F = 



Finally, we arrive at 



d 2 ,\ 
+ (-" 



E + m cos A 



)- 



dr 2 r 2e + (E + mcos A)r dr 

9 9 2eE e 2 — v 2 2em sin A + v cos A(E + mcos A) 

+E 2 — m -\ 1 = ; — 

r r z r [2e + (E + m cos A)r\ 

Let us introduce special designation for the additional singular point 

2e 



F = 



E + m cos A 



= R 



then we obtain 



£_ 1_ 1 d 

dr 2 r r — R dr 



o 9 2e£ e — v mR sin A - v cos A 
+£ 2 - m 2 + + =— + - 



r(r-R) 



F = 



(96) 
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After changing the variable x = r/R, it reads 



d 2 F 
dx 2 



+ 



X x — 1 



dF 

dx 



+ 



u 2 - P 2 



— v cos A + mi? sin A 2eRE — mR sin A + ^ cos A 



F = 0. 



x — 1 x 
Let us search solutions in the form F = x a e bx <p(x), the function ip obeys 



dx 2 



+ 



2a + 1 



x 



+ 26 



1 



x- 1 



dip 
dx 



+ 



6 2 + (E - m 2 ) R 2 + = 



J - ' - 2 - 2 '■■!)■ //cos .1 mi? sin .1 



x- 1 



a + b + 2ab + R(2eE - m sin A) + i/ cos A 

x 



V? = 0. 



When a, 6 taken according to (below we will use underlines values) 



a = +y/ 



v 2 — e 2 . 



■\Jv 2 — e 2 



b = +yjm 2 — E 2 R, - \J m 2 - E 2 R 
the above equation becomes simpler 



dx 2 



+ 



26 + 



2a + 1 



x 



x- 1 



dx 



+ 



a + 6 + 2a6 + 2eRE — mR sin A + v cos A a + 6 + v cos A — mi? sin A 



x 



x- 1 



In can be recognized as a confluent Heun equation for G(a, (3, 7, S, 77, z) 

G. + L + i±£ + l±2V 



+ 



z z — 1 

la + a/3-/3-/37-7-2r/ 1 a + cry + /? + /?7 + 7 + 2<5 + 277 



+ 



z-1 



G = 



with parameters 



a = 26, j3 = 2a, 7 = -2 , 
5 = 2 e£i? , 77 = 1 + mi? sin A - 2 eER - v cos A . 
Let us use the known condition for polynomial solutions 

/3 + 7 + 2^ 



n + 



it results the energy quantization rule 



a , n = 0, 1, 2, 



a = +\v 2 



b = - Vm 2 - E 2 R, 



eER = (n + y V - e 2 ) V ' m 2 - E 2 R , 



so we arrive at 



E = 



m 



l + e 2 /(n + vV 2 -e 2 ) 2 



which coincides with the known formula for energy levels. 
Let us consider the case 2). Eq. (5) takes the form 



— H — cos A — m sin A ) F = [ — 

dr r 



v sin A + e 



- 2m cos A) G , 



d v . . , \ e - usmA ^ 

cos A + m sin A ) G = F . 

dr r ) r 



One can obtain a second order equation for G(r): 



, d v . . . d v . 

f— H — cos ^4 — m smi r cos ^4 + m smA 

dr r dr r 



+ (e — usmA)( h 2m cos A) 



G = 0, 



from whence it follows 



d v 



( d 2 



cos A + m sin A + r ( — - + 

ar r 



^ cos A , v 



+ (e — ^sinA) 



\dr 2 
e + z/ sin A 



(- cos ^4 — m sinA) z 

r 



+ 2m cos A 



G = , 



or 



d} 



1 d 



-r^r H r ~ m sin ^ + 

ar z r dr 



e 2 — v 2 2me cos A m sin ^4 



+ 



+ 



G = 



And further, taking into account identity cos A = E/m, we arrive at 



d 2 1 d o 

t^t + --r- + -E — m + = — + 



e — v 



2eE Vm 2 - E 2 



+ 



dr 2 r dr 

Making the change of variables x = 2 v 7 rn 2 — E 2 r, we get 



G = 



d 2 G \ dG ( 1 v 2 - e 2 1 m 2 - E 2 + 2£e Vm 2 - -E 2 



<ix 2 x dx 



+ -T- 



c 2 + 2 



(m 2 — E 2 ) x 



G 



Let G(x) = x a e (p(x), the function 99 satisfies 



x-r-^ + (2a + 1 + 26x) -f- 
ax z ax 



+ 



(b 2 



)X + 



a 2 — v 2 + e 2 



+ 2ab + b + 



lm 2 -£ 2 + 2£e Vm 2 - E 2 



x 



rn 



2 — E 2 



<P = 0. 



When 



a = \J v 2 — e 2 , 



2 ' 



we get 



Ee 



d 2 (p dip 

x -r-TT + (2a + 1 — x) — a 

dx 2 K ' dx \ Vm 2 - E 2 



ip = 



which is a confluent hypergeometric equation 
x Y" + (7 - x)Y' - aY = , 



a = a — 



7 = 2a + 1 . 



\Jm 2 — E 2 ' 

Solutions become polynomials if a = — n, n = 0,1,2,..., this provides us with the energy 
spectrum 



E = 



m 



^l + e 2 /(n+^ 2 ^) 2 
In turn, from (12) it follows a second order equation for F{r) 



(13c) 



d v . . . 

— cos A + m sin A 

dr r j v sin A + e + 2m cos ir Idr r 



- — | — cos A — m sin A I F 



v sin A — e 



that is 



, d v . d f . ... 

(- cosA + m sin ^4) (- — I — cos ^4 — m sin ^4) 

dr r dr r 



+ I — In 



dr i/ sin A + e + 2m cos Ar J dr r 
e + z/ sin A + 2m cos 4r e-v sin ^4 



, d v . 
(— H — cos A — m sin ^4) 



After simple transformation, we arrive at 



d 2 n 
+ 



2m cos ^4 



F = . 



d 



dr 2 \ r sin ^4 + e + 2m cos Ar J dr 
2m cos ^4 



f sin A + e + 2m cos A 



-(- 

r Vr 



cos A — m smA , 



+£ 2 - m 2 + 



z/ 2 2eE — msin^ 



+ 



F = 



With special notation for the additional singular point 



(14a) 



9 



R = - 



e + v sin A 
2m cos A ' 



it reads shorter 



dr 2 



+ 



r r — R I dr r — R 



d m sin A v cos A , 1 1 . 
+ — (— —--) 



R y r-R r 



+E 2 - m 2 + 



v 2 2eE — msmA 



+ 



F = 



In the variable x = r/R, it looks simpler 



rf 2 F 
dx 2 



+ 



X x — 1 



dx 



+ 



zv 2 -e 2 



— v cos j4 + mi? sin A R(2eE — m sin A) + v cos A 



x-l 



F = 0. 



Let F = x a e bx <p(x), the function 99 satisfies 



dx 2 



+ 



2a + 1 oI 1 
+ 26- 



x- 1 



d</7 

dx 



+ 



2 /j-,2 2\ D 2 a — 1/ + e a + 6 + v cos ^4 — raR sin ^4 



6^ + (E 2 - m 2 ) R 2 + 



When a and 6 are 



it reads 



x — 1 



a + b + 2a6 + J?(2e.E - m sin A) + v cos ^ 



x 



V? = 0. 



a = +\/ zv 2 — e 2 , —\f v 2 — e 2 , 
b = +Vm T ^E 2 R, -Vm 2 -E 2 R, 



d^P + ( 2b+ 2a + 1 
dx 2 V x 



1 



d<£> 



x — 1 / dx 



+ 



a + 6 + 2ab + i?(2e75 - m sin A) + 1/ cos A 



a + b + v cos A — mR sin A 
x-l 

which is a confluent Heun equation for G(a, f3, 7, 5, 77, z) 



t> = 0, 



Z Z — 1 



G' 



+ 



1 a + a/3 - f3 - f3j - 7 - 2r/ 1 a + cry + /? + £7 + 7 + 25 + 2^ 

2 z + 2 z"^ J 



G 



10 



(15c) 



with parameters 



a = 26, (3 = 2a, 7 = -2, 
5 = 2eER, r] = l + mRsmA-2eER-iscosA. (15d) 
Imposing the known condition for polynomial solutions 

5 = - (n+ P + 1 + 2 ] a> n = 0,1,2,... (16a) 



2 

we derive the energy quantization rule 



+ yV 2 - e 2 , b = -y / m 2 -E 2 R, 



eER = (n + \Av 2 - e 2 ) \/™ 2 - E 2 R 



m 



E = , , (166) 

/ l + e 2 /(n+ V^-e 2 ) 2 



which coincides with the known result. 

It should be emphasized that confluent Heun equations in cases 1) and 2) formally coincide, 
however all parameters are different in fact: 

1) a = 26, (3 = 2a, 7 = -2, 

5 = 2eER, j] = 1 + mR sin A- 2 eER- v cos A, (17a) 

2) a = 26, /3 = 2a, 7 = -2, 

5 = 2eER, r]=l + mRsmA-2eER-iscosA, (176) 

where 



2e e e 2 

R = -— , smA=-, cosA = \ l-^; (18a) 

E + m cos ^4 v V v z 

e + usinA t E t I E 2 

R = , cosA= — , smA = \l = . (186) 

2E m V m 2 y 1 

Let us turn again to radial equations in presence of Coulomb potential 

(-f + - )/ + (E+- + m)g = 0, 

ar r r 

(j- - -)g - {E + - - m)f = 0. (19) 

ar r r 



Excluding the function g, one gets 



d 2 f e df 
— - H + 

dr 2 r(Er + e + mr) dr 



e(e 2 - v 2 ) 
r 2 (Er + e + mr) 



11 



e (E + m) (3E — m) r(E — m)(E + m 



Er + e + mr 



+ 



Er + e + mr 



f = o. 



After changing the variable 

eq. (20) takes the form 

x d 2 f 1 df 
dx 2 x — 1 dx 

Separating two factors 



x = — - 



(E + m)r 



+ 



e 2 (Ex — mx - 2E) e 2 — v 
E + m ' x 

f(x) = x A e Cx F(x) , 



x- 1 



f = o. 



for F one derives 



dx z \ x x — 1 J dx 



C 



2 e 2 (E - m) 



E + m 



A 2 + e 2 - v 2 A + C + 2^C - 2Ee 2 /(E + m) + u A + C + v 



x- 1 



When A, C are taken as (bound states are of the prime interest) 



E + m 
A 2 + e 2 - v 2 = 



C = +e 



m — E 
m + E ' 



eq. (23) becomes simpler 



A = +v / z^ 2 - e 2 



+ 



dx 2 V x x — 1 J dx 

A + C + V + 2AC -2Ee 2 /(E + m) A + C + v 



F = Q, 



x x — 1 

which is the confluent Heun equation for F(a, /3, 7, 5, rj; x) 

f_ F+ ( a + P±± + l±±) 

dz 2 \ z z — 1 / dz 

( 1 a + aj + (3 + + 7 + 25 + 2r] l a/3 + a - /3 7 - /? - 7 - 2tA 
+ \2 z-l + 2 z ) 



with parameters determined by 



12 



a - 2( ' - +2, a/— , p = 2A = +2 y/v 2 - e 2 , 



7 = -2, 6 



m + E ' 

2£e 2 



E + m 

The known condition to reach polynomials is 



rj = 1 — v + 



2£e 2 
£" + m 



<5 = —a ( n + 



7 + /3 + 2 



It results in 



2Ee 2 



-2e 



m — E 



or 



from whence it follows 



E + m \ m + E 

Ee 



( n + \J v 2 — e 2 ) , 



yjm 2 — E 2 



= N , N = n + ^u 2 - e 2 ; 



m 



VI + e 2 /iV 2 ' 

It is the exact energy spectrum for hydrogen atom in the Dirac theory. 



(266) 
(27a) 

(276) 
(27c) 

(28) 



Supplement A. Coulomb problem in spaces of constant curvature 

In the Lobachevsky space H3, the problem of Dirac particle in Coulomb field reduces to a radial 
system (let us specify the case 5=1) 



d v 
+ 



/ + (E + 



( 



dp sinh /3 tanh P 

d v , _ e 

) 9 ~ (E + 



"dp sinh P tanh P 

From this it follows a second order equation for / 

d 2 f 



+ 



+ m) g = , 
- m) f = . 

df 



{Ala) 



+ 



{E + 



dp 2 sinh P [e cosh P + (E + m) sinh /3] dx 
2 



tanh /3 



o 2 1/^ + 1/ cosh /3 
! ^ m — h 



sinh 2 P sinh 2 P e cosh P + (E + m) sinh /3 



/ = o. 

(A.lb) 



Let us apply a linear transformation with a(P)b(P) — c(P)d(P) = 1: 

/(/3) = a F(/3) + c G(/3) , g{p) = d F(P) + 6 G(/3) , 
F(/3) = 6 /(/3) - c g(P) , G(P) = -d f(P) + a 5 (/3) . 
After combining eqs. (A. la) in the same way as in the flat space model we obtain 



(A.2) 



13 



d , / / , v(ba + cd) . „ e . , , . „ e 

Tr ha+ " d+ ^^ +{E+ ^ +m)hd+{E+ ^rr m)ca . 



b'c - be' 



sinh (3 



26c - (£ + 



tanh (3 

e 



+ m)b 2 -{E + 



tanh (3 



— m)c 2 



d u(dc + ab) , „ 



+ m) fed - (£ + 



tanh /3 



-d'a + da' + . - 2ad + {E + — 

sinh p tanh p 



+ m)d 2 + (E + 



tanh /3 



— mm 



G, 

G 
F . 



When the above transformation is orthogonal 
eqs. (A. 3) become simpler 



a c 




cos A/2 


sin A/2 


d b 




-sin A/2 


cos A/2 



d z^cosA . , , 

— + . - -m sin A ] F 
dp sinh p 

d vcosA . . 

+ m sin ^4 G 



d/3 sinh /3 



A e cosh /3 + i/ sin ^4 
~ Y sinh (3 

A' e cosh f3 — v sin ^4 
+ T + sinh (3 



— E — m cos A j G , 
+ E — m cos A] F . 



Let ut translate equations (Ala) to a new variable 



eqs. (Ala) will take the form 



£/+-/- 

az z 

d ^ 

-1-9- -9 + 

dz z 



, 

tanh — = z ; 
2 



2(£ + m) , e(z 2 + l) -, 

Z 2 -l + z(z 2 -l) , y 

2{E-m) e(z 2 + l) 



z 2 _ ! 



+ 



z(z 2 - 1) 



/ = 0. 



or differently 



dj_^_ V j_^fe_^_— E — e — m ^ E — e + m\ ^ 
dz z \z z — 1 z + ly 

d f 6 + ~^ ~ e + m E — e — m\ ^ q 
dz z \z z — 1 z+ly 

From (A. 6) it follows a second order differential equation for / 



d 2 / 

dz 2 



+ 



1 1 1 

- + T + 



ez + E + m 



df 

z ' z — 1 ' z + 1 " ez 2 + 2(E + m) z + ej dz 
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+ 



2Ee 2 - (E + m) v -(E + e) 2 + m 2 + v (E - e) 2 - m 2 - v e 2 - v 2 

2 + " + " + q 



ez 



z-1 



z + 1 



{E + e) 2 -m 2 (E-ef-m 2 2 v [ez{E + m) + 2 (E + m) 2 - e 2 } 1 



(z-1) 2 (z + 1) 2 e[ez 2 + 2(£ + m)z + e] 

It is convenient to use shortening notation 

E + m 



f = 



= a . 



(A.7a) 



Eq. (A 7a) has 6 singular points 



, oo , ±1, 
zi, 2 = -a ± Vo" 2 - 1 , (z x z 2 = 1 , zi + z 2 = -2cr), 



and it reads 



d 2 / 

dz 2 



+ 



1111 1 

- + 7 + 



z Z — 1 Z + 1 Z — Z\ z — Z<l 



dz 



+ 



AEe-2av (E + e) 2 - m 2 - v (E - e) 2 - m 2 - v e 2 - u 2 



z - 1 



z+1 



(E + e) 2 -m 2 (E — e) 2 - m 2 A B 



(z-l) 2 (z + lf 



Z — Z\ Z — Z2 



/ = 0, 



where 



a z + 2a 2 -1 A B 
2v- — r = + 



(Z - Zl) (Z — Z 2 ) Z — Zl Z — Z 2 

^ = 2 ^z 1+ 2. 2 -l ^ = 2 ^ 2 + 2. 2 -l 



zi - z 2 



Z2 ~ Zl 



Let us introduce notation 

C = (E + e) 2 - m 2 , D = (E - ef - m 2 , AEe = C - D , 
then eq. (A.7b) can be presented shorter 

d2 f + f i + _J_ + 1 1 1 # 



<iz 2 \ z z — 1 z + 1 z — zi z — Z2 J dz 

(C-D-2ov C-v D-u e 2 -v 2 

+ T + ^T + — 

\ z z — 1 z + 1 z^ 

C D A B \ £ n 

(z-1) 2 (z + 1) 2 z-zi z-z 2 ) 
Behavior of its solutions near 3 singular points 0, +1, —1 is given by 



(A.7b) 



(A.8a) 



(A.8b) 



(A.9) 
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x = +1 , /" + -/' + Cf = 0, f~z a ,a = ±V=C ; 

x — 1 

x = -l, f"+-^—f' + Df = 0, f~z? , (3 = ±V=D; 
x + 1 

x = 0, f" + lf' + (e 2 -^)f = 0, f~z M , M = ±s/v*-# 



Besides, near the singular point z = oo we have 



x = oo 



XX X* 

remembering that A + B — 2au = we get 

a~/, N = ±y / C + D + e 2 - u 2 . 
Let us search solutions in the form 

M( 

J = * 

eq. (A. 9) leads ro 



f = x M (z-l) a (z + l)e<p; 



dz 2 



+ 



2M + 1 2a + 1 2/3 + 1 



+ 



Z — 1 Z+1 Z — Zl Z — Z 2 

M 2 + e 2 - z, 2 a 2 + C /3 2 + D 

n + 7 TT^ + 



dz 



(z-1) 2 (z + 1) 2 
C-D-(a-f3) (2M + 1) - 2<r(v + M) 



M + a/2 + /3/2 - C + v + 2Ma + af3 
M + a/2 + /3/2 - D + v + 2M/3 + a/3 



+- 



z — Zl 



a 



/3 



z+1 

1 



Z\ — 1 Zl + 1 Zl y Z - Z 2 



a 



/3 MY 



z 2 — 1 z 2 + 1 z 2 / 



Requiring 

M 2 = ±yV - e 2 , a = ±v /Z C, /3 = ±^/^D , 
we arrive at a differential equation with six singular points (0, oo, +1, — 1, zi, z 2 ) 



dz 2 



+ 



2M + 1 2a + 1 2/3 + 1 



z — 1 z+1 Z — Zl z — z 2 



dz 



(AlOa) 



(A106) 



= . 
(Alia) 
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+ 



C-D-(a-P) (2M + 1) - 2a {y + M) 



M + a/2 + p/2 - C + v + 2Ma + a/3 
z-l 

M + a/2 + (3/2 - D + v + 2M/3 + a/3 



z — Z\ 



A 



z + l 

a 



/3 M\ 



+ 



z - z 2 



B - 



Z\ — \ Z\ + 1 Zi / 

a f5 M 



LP = . 



Z 2 — 1 Z 2 + 1 Z2 . 

(A116) 

Let us again try to perform a linear transformation of the type (A. 2) - which results in 



b'a + c'd + - (ba + cd) + [ - + 

dz z 



-E — e — m E — e + ra 



z-l 



+ 



z+l 



bd 



e — E — e + m E — e — m 

+ z-i + ^n- lcn 



b'c -be' - - 26c 

z 



e —E — e — m E — e + m\, 

~ + 1 + Zl — ' b 

z z — I z+l 



e —E — e + m E — e — m 

_ + 1 — + — in — 

z z — 1 z + l 



G 



dz 



+ d'e - a'b-- (dc + ab) - ( - + 

z \z 



e —E — e — m^E — e + m 



z - 1 

E — e + m E — e — m 



z+l 



bd 



e 

- + 

z z — 1 



+ 



z+l 



ca 



G 



-d'a + da' + - 2ad + ( - + 



-E — e — m E — e + m 



z 



z - 1 



+ 



z+l 



, e — E — e + m E — e — m 

+ ( - + : + — j— 

. z z — 1 z+l 



When the transformation is orthogonal, eqs. (A. 12) become simpler 



(A12) 



d v 



+ - cos A - - + 



dz z 



e —E — e — m E — e + m\ sin^4 



z - 1 



+ 



z + l J 2 
41 

F 



e —E — e + m E — e — m\ sin A 
H z + z-l + ^+l-J I" 



1 i> . . / e — £ — e — m E — e + m 
sin A- - H 1 

I z \z z — 1 z+l 



A 



cos 



e — — e + m E — e — m 

- + 1 + ZT— 

z z — 1 z + l 



■ 2^ 

sin - 



G, 
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d v f e —E—e—m E — e + m\smA 

dz z \z z — 1 z + 1 / 2 



e —E — e + m E — e — m\ sin A 
~ I + z - 1 + z + l J 2 



G 



+ - 



A' 



v 



sin ^4 



2 



2; 



+ 



-E — e — m E — e + m \ 

H : I sm z 



z 



1 



e — E 1 — e + m E — e — m 

+ ( - + ; + —, — 

z z — 1 z + l 



z + l 

COS - 



A 
~2 



After simple algebraic manipulation we arrive at the system 



d vcosA msinA msinA 
dz z z — 1 z + l' 



A' e + v sin A 



-E — e — m cos ^4 I? — e + m cos A 



z z — 1 z + l 

d vcosA msin^ msinA N 
dz z z — 1 z+l 



G, 



G 



A' e — vsuiA —E — e + mcosA E — e — m cos A 

+ + ; + ; 

2 z z-1 z+l 



F . 



In the limit of flat space model, z << 1, they lead to the known formulas (1.5). 
Let us try the variant (see Section 1.10) 



in this case eqs. (A. 13) read 



d_ 

dz 



sin A = ejv ; 



v cos A m sin A m sin A 



2v 



■ sin A — 



z z — 1 z + l 

-E — e — m cos A E — e + m cos A 



z z — 1 

v cos A m sin A m sin A 



z+l 



G, 



d_ 

dz 



Excluding the variable F, we get 



+ 



z + l 



(2mcosA-2£) - 2ez 
(*-!)(* + !) 



(A13) 



(A14a) 



(A146) 



d z^cosvl msinA msinA\ (z — l)(z + 1) 
dz z z-1 z + l J (2m cos A - 2E) - 2ez 

( d vcosA msinA msinA\ . 

x 1 ) G 

\dz z z — 1 z+l 



2u . . —E — e — m cos A E — e + mcosA\ 
sin A — G , 

z z — 1 z+l 
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(A.Uc) 



from whence it follows 



^ (d^ (z-l)(z + l) 
dz 2 V dz (2m cos A — 2E) - 2ez 



d z^cos^ to sin A to sin A" ^ , 
dz z z — 1 z + l 



/ 1/ to sin A to sin ^4 

- -cos ,4 + — - 

\z z — 1 Z+1 

2tocos,4 - 2£ - 2ez /2i/ 
(z-l)(z + l) \z 



G — 



v cos A to sin A to sin A 



+ 



: z — 1 z+1 

—E — e — to cos A E — e + to cos A 



z-1 



2 + 1 



G = 0. 



It is convenient to introduce a special designation for the singular point zq = (mcosi — E)/e, 
then 



+ 



+ 



1 



1 



+ 



d 



z-1 z+1 



z + 1 z — zq J dz 
v cos A to sin A to sin A 



z-1 



+ 



z+1 



G 



— ( - cos A 

z 



m sin A to sin A 



—e 



1 - -sp 
z - 1 



+ 



z 

z + 1 



1 



z+1 



G + 



v cos A to sin A to sin A 



+ 



(z-1) 2 (z+1) 2 



2v -E 

— sin^H 

z z — 1 



e — to cos A E — e + to cos ^4 



z + 1 



With the use of notation 

to cos A = c , to sin ^4 = s , ^ cos ^4 — > v , 
and taking in mind identity usuvA = e, eq. (A.14d) is written shorter 



G 

G = 0. 

(A.Ud) 

(A15) 



+ 



<iz 2 ~^ \ : 

1 

+ 



1 1 

+ 



1 



z — 1 z + 1 z — Zo 

2 

+ 



d_ 

l'z + l z — zo J dz 

1 \ I V 

z 



+ 



V s s 
z z — 1 z + 1 



+ 



z-1 z+1 

s s 



1 - z 1 + z 
-e I + 



z-1 



z+1 



2e 



+ 



z 2 (z-1) 2 (z+1) 2 
— E — e — c E — e + c 



z-1 



+ 



z+1 



G = 



(A16) 



Equation (A. 16) reduces to 



1 1 

+ 



1 



d 



^G + 

dz 2 V z — 1 ' z + 1 z — zq J dz 
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+ 



1 —v + Avszq - 4z 2 e 2 1 ZqV + 2z s - v v{y — \) 



+ 



Z- Z Zq(1 + Zq){zq - 1) 



1 



z + 1 

1 



{+v — s 2 — 2vs + e 2 + 2e 2 zo — czqc — ez$E + 



z-1 



{—v + s 2 — 2us — e 2 + 2e 2 zo + ezoc + ezoE 1 + ) 

1 — Zq 



-s 2 + (1 - z ) (e 2 + ec + eE) -s 2 + (1 + z ) (e 2 - ec - e£) 

H 7 TTo "I - 



G = 0. 



Let us specify behavior of G(z) near the singular points 

0, oo, +1, -1 . 



Near the point z = we have 

d 2 



dz 2 zo dz 
therefore G behaves in accordance with 



G ~ z a , a(a - l)z a ~ 2 + — + a z a_2 z - 1/(1/ - I)* 



1 



a-2. 



_a-2 
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a(a — 1) = z/(z/ — 1) , a = +zv, 1 — v ; 



positive values a = +zv correspond to bound states. 
Near the point z = 1, we have 



that is 



^ 2 C| 1 ^ | - g 2 + (l-2o)(e 2 + ec + e^) C _ Q 
a!z 2 z — 1 dz (z — l) 2 



G~(z-1) M , M = ± V* 2 " (1 " 20) (e 2 + ec + e£) 



because z = tanh (/3/2), the sign minus in (A. 186) corresponds to the bound states. 
In the same manner, consider the vicinity of the point z = — 1: 

^ 2 c , 1 d | -s 2 + (1 + go) (e 2 - ec - eE) Q _ Q ^ 



<iz 2 z + ldz 



{z + lf 



that is 



G~(z + l) JV , iV = iyV - (1 + z ) (e 2 - ec - eE) ; 



(A17) 



(A.18a) 



(A.m) 



(AA8c) 



because z = tanh (/3/2), the singular point z = — 1 is not physical one. 

For brevity, let us introduce special notation for coefficients at the for singular terms 



K 

z- Zq 



Ko 

z 



x — 1 ' X + 1 ' 
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then eq. (A. 17) reads shorter 

dz 2 



G + 



d 



z — 1 z + 1 z — zq J dz 



1 1 

+ 



G 



+ 



K K v(v - 1) K+ K_ 

— + ^ „ ' + ±7 + 



M 2 



TV 2 



z - z 



z+l z-l (z-l) 2 (z + 1) 2 _ 



G = 



Let us introduce the following substitution 

G = z a (z-l) M (z + l) J >; 
after simple calculation we arrive at 



dz 2 



+ 



+ 



n 2 -iV 2 
'(z + 1) 2 



2a 2m + 1 2ra + 1 1 

z Z — 1 Z + l Z — Zo 

a (a - 1) - v(y - 1) m 2 — M 2 

m 



+ 



+ 



z - z 



dip 
dz 



n 



z z - 1 z + 1 



1 / ^ a[-2mz + 1 + 2nz ] \ 

1 / T , 1 , x , 1 m-fzn + 1) 

+ — ( A - + 2 (2o + " ) < 2 '" +1) + 2 4^T i 



ip = 0. 



Requiring 



we arrive at 



a(a - 1) = v(y - 1) , m 2 = M 2 , n 2 = N 2 , 



dz 2 



+ 



2a 2m + 1 2n + 1 1 
— + r + ^^ 

z z — 1 Z + l Z — Zq 



dz 



+ 



z - z 



a m n \ 1 , , 

# r — T + - ( ^ 4 

z z - 1 z + 1 J z 



a(—2mzQ + 1 + 2nzo) x 



1 



z-l 
1 

"z + l 



+ - (2a + n) (2m + 1) + - K ' 
2 2 Zo — 1 

1 , , . 1 n(zo — 1) 

-2(2a4m)(2n + l)--^_i 



V? = 0. 



In the Riemann space, the Kepler problem can be treated in a similar way. 
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